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ABSTRACT 


A study of heat transfer rate of spherical copper and alumina nanoparticles in water and ethylene glycol based fluids was carried out. 
A modified thermal conductivity model in conjunction with steady state momentum and energy equations in spherical coordinates 
were put into dimensionless form and solutions used to determine the skin friction, heat transfer coefficient and thermal 
conductivity as well as viscosity. The modified model incorporates Brownian motion and varied sphericity to observe the effect of 
temperature and other material parameters on the velocity and temperature profiles of the fluid. Using numerical values, it was 
shown that the nanoparticle volume fractions, the diameter and Prandtl number, not only enhanced the thermal conductivity of 
nanofluids but also the velocity and temperature profiles. It was also observed that the Brownian motion which is temperature 
dependent was actually a weak factor in enhancement of thermal conductivity. The effect of other parameters as well as calculation 
of mass flux and mean temperature was determined. 
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1. INTRODUCTION 


Generally, nanofluids are formed when some quantity of nanoparticles and base fluids are mixed. The nanoparticles has an average 


diameter between 1nm to 100nm.The nanoparticles are mainly oxides, metals and carbides while the common base fluids are water, 
ethylene glycol and other traditional heat transfer fluids. Nanofluids exhibit enhanced properties when compared to conventional 
heat transfer fluids and other metallic nanoparticles. The larger relative surface area of nanoparticles, compared to those of 
conventional particles, appears to justify significantly the enhanced heat transfer capabilities and also should improve the stability of 
the suspensions. Also, nanofluids enhanced abrasion-related properties as compared to the conventional solid/liquid mixtures. 
Successful use of nanofluids will support the component miniaturization by enabling the design and production of smaller and 
lighter heat exchanger systems. In his pioneer work, Choi [1], described the mixture as nanofluid with increased thermal conductivity 
of the base fluids and their convective heat transfer rate. The mixture or nanofluid when compared to the base fluid, it is observed 
that changes occur in viscosity, thermal conductivity and density. Nanofluids can therefore be described as new class of enhanced 
heat transfer fluids formed as a result of dispersion of nanometer-size particles in base fluids. Of the basic physical properties of 
nanofluids, thermal conductivity is the most important due to its application in heat transfer. The enhanced heat transfer was 
established because; nanofluid thermal conductivity is a function of nanoparticle size, shape, volume fraction and spatial distribution. 
Others are base fluid type, temperature and pH value (Yu et al [2], Eapen et al [3]). Nanofluids uses are abounded in various 
electronic equipment, energy production, power generation and air conditioning and production. Serious study in convective heat 
transfer using suspensions of nanometer-size particles in a chosen base fluid started only over the past decade. The study is either 
experimental or theoretical and a combination of both. Keblinski et al.[4], in a study, outline and examined the properties of 
nanofluids and future challenges. In spite of the several successes made, the research and development of nanofluids is seen as 
rapid but still hindered by several factors such as the lack of agreement between experimental and theoretical results, poor 
characterization of suspensions, and the lack of theoretical understanding of the mechanisms. This trend rather increased the study 
of naofluids and its applications. According to some proposed classical nanofluid models (Xie et al [5], Hamilton and Crosser [6], 
Jeffrey [7], Davis [8], Wang et al [9], Koo et al [10], Jang and Choi [11]), effective thermal conductivity of mixtures can be calculated 
and other parameters of the system also determined. However, results of some researches (Patel et al [12], Das et al [13], Xuan et al 
[14], Kumar et al [15], Bhattacharya et al [16], Putnam et al [17], Koo and Kleinstreuer [18 and 19]) strongly argued that the non 
inclusion of Brownian motion effect which is temperature dependent to the earlier stated models is a minus and therefore not a 
holistic approach to the determination of thermal conductivity enhancement, hence their studies. Surprisingly, Wang et al [20], in 
their study, declared that the effective thermal conductivity enhancement due to Brownian motion is not necessary. The implication 
is that the effect of temperature on the enhancement of thermal conductivity is very minimal and can be discarded. Keblinski [21], 
opined that the heat transferred by nanoparticle diffusion contributes minimally to thermal conductivity enhancement. In all, errors 
arising from preparation of nanofluids include, heating and measurement processes and cleanliness of apparatus are some of 
experimental results differences in the various studies while over simplification and assumption of certain parameters are the causes 
of some classical models. As a result of these uncertainties, some of the results appear to be conflicting. The aim of this study is to 
include Brownian motion into our proposed model which is a modification of the classical models of nanofluids of Koo and 
Kleinstreuer ([18 and 19]). By this modification, different sphericities for different shapes nanoparticles in addition to other 
consideration which is an extension are made. 


Formalism 
A steady boundary layer flow of viscous, incompressible, Newtonian Spherical nanofluid is considered. Using the Boussineq’s 
approximation, the governing equations of the nanofluid in steady spherical coordinate is given as 


Hp { 2v, 4 Ov omy 
— “~+—~ |+ gf (1 -T, )=0 (1) 
Py r? r’ or’ or’? By ( Fl 
Ki 20f OT 
if 6) 6) _0 . 


(oC, : ror r Or’? 


Subject to the boundary conditions (Wang et al, [9]) 
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V,. ' = Vo 
r=R 
(3) 
V, r'=0 Vo 
r'=R = fy 
r (4) 
r'=0 0 


where V, is nanofluid velocity, V_ is characteristic velocity, r is radius of nanoparticles, T is temperature of nanofluid, I, is 


characteristic temperature, 7’, is ambient temperature and g is acceleration due to gravity. 


In this work, according to Hamilton and Crosser model [6], effective dynamic viscosity ratio and effective thermal conductivity ratio 
which are valid for both spherical and non spherical shapes nanoparticles are stated as 


Pw = (1+a9+b¢?) (6) 


iy 
ky — k, +(n-Ik, +(n-1k, —k, 
k,  k, +(n—1k, —{k, -k, 


(6) 


3 
The n is the empirical shape factor given by m1 = — where y is the sphericity. 


According to the work of Tiwari and Das [22] and Asma et al [23], density of nanofluid (Pr ), thermal expansion due to 


temperature of nanofluid (Buy ), specific heat at constant pressure of nanofluid are respectively 


Pru =(1-¢)p, +p, 
By =(1-9)B, + 9B, 7) 
(c, ), = (1 — oC, }, + dlc, ) 


v ys 
where @ is the nanoparticles volume fractions given by ¢ = ——~— = m—D>, m is the number of particles per unit volume 
Ve +tVv 
ef s 


and D is the average diameter of the particles, , and P, are the densities of the base fluid and solid nanoparticles, B; and 
B, are the thermal expansion due to temperature of base fluid and solid nanoparticles, and (c, ), and (c, i} are the specific 


heat at constant pressure due to base fluid and solid nanoparticles. a and b are constants that depend on the particle shape (Aaiza 


et al [24]). The thermo physical properties of alumina (Al,O;) and copper (Cu) nanoparticles as well as water (H,0) and 


ethylene glycol (C,H .O,) as base fluids are presented in table 3. 


Table 1: Sphericity Y and empirical shape factor for different shapes nanoparticles (Aaiza et al [24]) 


Model Platelet Blade Cylinder Brick 
y 0.52 0.36 0.62 0.81 


n 5.76923 8.33333 4.83871 3.70370 
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Table 2: Constants a and b empirical shape factors (Aaiza et al [24]) 


Model Platelet Blade Cylinder Brick 
a 37.1 14.6 13.5 1.9 
b 612.6 123.3 904.4 471.4 


Table 3: Thermo physical properties of Al,O, and Cu nanoparticles, C;H,O, and HO (Aaiza et al [24]) 


Prepery H,O C,H,0O, Al,O; Cu 

p (kg/m *) 997.1 1115) 3970 8933 
k, y (W/mk) 0.613 0.149 40 401 

Lag (m?s"') 0.00089 0.001095 0.0000029014 0.00046 
B xl0°(k) 21 6.5 0.85 1.67 


The Hamilton and Crosser model [6], shows that the effective thermal conductivity of nanofluids depends on the thermal 
conductivity of the spherical solid, the base fluid and the volume fraction of the solid parts. However, to effectively determine the 
thermal conductivity of nanofluids, Patel et al [12], observe that, the Brownian motion effect on nanoparticles at the molecular and 
nano-scale levels may be a key mechanism governing the thermal behaviour of nanofluids. In a similar experimental study 
conducted by Das et al [13], they find out that thermal conductivity of nanofluids depends strongly on temperature and that this fact 
should be considered in theoretical models. The Hamilton and Crosser model ignored Brownian motion which is temperature 
dependent; it is therefore defective in determining the effective thermal conductivity of nanofluids. In a related development, studies 
that incorporate Brownian motion were reported. Few of the works include Xuan et al [14], which included the Brownian motion into 
the model proposed but was criticized of being too weak in temperature dependence and also at variance with experimental data of 
Das et al [13]. Also, Kumar et al [15], proposed a model to account for the large enhancement of thermal conductivity of nanofluids 
and its strong temperature dependence, but its use showed that it is not suitable for high concentration of particles. Another study 
carried out by Bhattacharya et al [16], developed a technique to determine the effective conductivity of a nanofluid using Brownian 
motion simulation. Although the model showed good agreement with other works, its sphericity is just 1. The same issue of 
sphericity being 1 was the weak point of the model proposed by Putnam et al [17] which is a modification of Maxwell's model. 
Above all, Koo and Kleinstrever ([18 and 19]), proposed a model of nanofluids which among other factors included the effect of 
particle size, particle volume fraction, temperature and properties of the base fluid. The model is stated as 


Ky Ky + 2k, + lk, =k +5 0°eleC,), | 1.9) (8) 


k,  k,+2k,—(k,-k,}p 


where & and hs are respectively, related to particle motion and Stefan- Boltzmann constant. The function f(T.¢) can vary 


continuously with particle volume fraction as f(T, ¢) = (- 6.04¢@ + 0.4705 WV + (1722 30¢ = 134.63). Comparing equation 
(6) with equation (8), it is realized that the first part of equation (8) is a special case with sphericity 1 of equation (6). To tackle the 
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effect of heat transfer, sphericity and temperature dependent Brownian motion of Cu and Al,O, nanofluids, equation (8) is 


modified to take the form 


ky _k,+(n-1k, + (nk, -k, 
k,  — -k, +(n—1k, —{k, -k, 


Ss 


+5x10*€¢(pC, ) 


= f(T.9) (9) 


2. DIMENSIONAL ANALYSIS 
To effectively tackle the governing fluid flow equations, dimensional homogeneity of the governing equations using the Buckinham- 
7 — theorem is stated 


pS see a — 8B, (T-T ey 


D D v,p,D) T-T,’ ky ve 


Also, substitute equation (9) into equation (2) and equation (5) into equation (1), the resulting equations transformed into 


(+agrbg?)(2u 40 4) (a gle geo a0 
Re r- r Or Or’ B, ’ 


6(2060 00 
+ =0 (11) 
Pri r ér or’ 


Subject to the boundary conditions 


u| f= 0 
7 (12) 
u| = 0 
=1 
r=l (13) 
r=0 =0 


where Re is Reynolds’ number, Pr is Prandtl number, G, is Grashofs number, @ is dimensionless temperature, u is dimensionless 


velocity and r is dimensionless radius of the nanoparticles. 


Ms suto “Eglo, ). 


_ k,+(n-1k; +(n-Dlk, - 


k,+(n—Dk, -(k, “oy 


Simplification of equation (11) results in 


f(T.) 


k 


C +5x10*E¢(pC,,) . (6.04 +0.4705)0!5 + (1722.39 - 134.63)0°° ] 


Ss 


200 06@ 
+ =(0) (14) 
Pr ror @r’ 


pase |. 7 


ARTICLE 


k, +(n-Ik, +(n—1\k, -k,.\b 
k, +(n—1k, —(k, —k, }b 


where @, Pr= 


Method of Solution 


Approximate g'° and @°° using Taylor's series expansion about 1 and neglect powers of 9 > 2 and simplify, equation (14) can 


be written as 


2 
(ate 2b-1a0;03)| — "2 e |a6 (15) 
ror Or 


where 


suo Eg(pC, ) “ (— 6.049 + oar) 


Ss 


G= 
2Pr 


suo‘ sec) rm9—13469) 


Ss 


a = 
: 2Pr 


The simplification of equation (15) and its solution as well as imposing the boundary conditions of equation (13), results in 
6° (3a, +a, )+ O(a, —a, +a,)—(a, + 2a, +2a, )r =0 (16) 


and finally reduced to 


(a, —a,+ a,)+ Ve, a ay y + 4(3a, +s \(@, 42a, + 2a,)r) 


vue 2(3a, + a,) 


(17) 


Since 9 > 0, +@ is chosen. 
To solve equation (10), equation (17) is put into it and the resulting expression is simplified as 


2 +4/B5 + 
Gare, B, B, + far _¢ 


r ror Or? PB; 


(18) 


The resulting solution of equation (18), after imposing the boundary conditions of equation (13) is given as 


pr Ware er) [alereay 1), aller) 


MO OB. ISB 158,82 2\ 158,82 


(19) 


where 
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(1-9)+9” |Rea, 


. a Troe B, =(a, ay +a,), B; = 4(3a, + a, ay + 2a, +2a,), Bs = 2(3a, +a,) 


The mass flux @ and the mean temperature 6, are obtained by evaluating the integral 


1 
oO = Jucryar (20) 
0 
and 
1 
[um ocr)ar 
0, = x [anes (21) 
[ucar 
0 


Heat Transfer Coefficient (Nu) 

It has been shown that the heat transfer performance is a better indicator than the effective thermal conductivity of nanofluids, 
particularly when such nanofluids are used as coolants and other functions in industries. This property necessitated the formation of 
nanofluid models to either complement the models of thermal conductivity or an improved determination of effective thermal 
conductivity. Following the models for determining the heat transfer coefficients of nanofluids (Polidori at al [25], Mansour et al 
[26]), the mean Nusselt number is calculated thus, 


60(r) = Nu= Bs (20) 
or r=0 2B, V B F 
Skin Friction 


The skin friction is given by 


éu(r) ee 2p; és 4{Vez +B.) | (21) 
ar |,» 38,B:(\B:} 158,8; 2 


3. RESULTS AND DISCUSSION 
Results 


Table 4: Numerical values of mean Nusselt number for various values of Prandtl number in H,O and C,H .O, based fluids 


when 


k, =1.380658 x10 JK" ,é =1,0 =25,d =15nm,¢= 0.06mm? 


Pr Nu(Al,O,)C,H,O, Nu (Al,O,)H,O Nu(Cu)C,H.O, Nu(Cu) H,O oO 
4 
0.71 — 3.60599 x107"” — 3.48063 x10” 2.07758 x107'* 2.06972 x107"” Sp 
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1.71 — 3.60597 x10~ — 3.48062 x10~” 2.07757 x10~? 2.06971 x10~”” 
2.71 — 3.60593 x10” — 3.48058 x10” 2.07756 x10” 2.06969 x107'? 
3.71 — 3.60588 x107” —~ 3.48053 x10~!” 2.07753 x10~? 2.06967 x10~? 


Table 5: Numerical values of mean Nusselt number for various values of nanoparticles volume fractions in H,O andC,H.O, 
based fluids when 


k, =1.380658 x10 * JK ',é =1,0 = 25,d =15nm, Pr = 0.71 


d(mm*) —-Nu(Al,O,)C,H,O, Nu(Al,O,)H,O Nu(Cu )C,H,O, Nu(Cu )H,O 


0.06 3.60599 x10~? 3.48063 x10~' 2.07758 x10? 2.06972 x10°? 
0.07 9.80488 x10" 9.41793 x10"? 5.6575 x10" 5.63316 x10 
0.08 165227 107? 1.57965 x10"? 9.54736 x10°'4 9.50156 x10-\4 
0.09 5.51207 x107"” 5.24614 x107 3.18945 x10” 3.17262 x10~"* 


Table 6: Numerical values of mean Nusselt number for various values of nanoparticles diameter in H,O andC,H,O, based 


fluids when 

k, =1.380658 x10” JK |, € =1,0 = 25,6 =0.06mm’, Pr = 0.71 

a(nm) Nu( Al,O,)C,H,.O, Nu(Al,O,) H,O Nu(Cu )C,H,O, Nu(Cu)H,O 
15 3.60599 x10 7” 3.48063 x10~? 2.07758 x107 2.06972 x10~'? 
17 4.08679 x10~? 3.94472 x10°- 2.35459 x107” 2.34568 x107” 
19 4.56759 x10” 4.4088 x10~”” 2.6316 x10” 2.62164 x10” 
21 5.04838 x107” 4.87289 x107? 2.90861 x107"” 2.8976 x10°” 


Table 7: Numerical values of mean Nusselt number for various values of nanoparticles Sphericity in HO and C,H .O, based 


fluids when k, = 1.380658 x10 JK 1, =1,0 = 25,¢ =0.06mm’, Pr = 0.71,d(nm) =15 


e Nu(Al,O;)C,H.O, Nu(Al,O,) H,O Nu(Cu )C,H ,O, Nu(Cu ) H,O 

0.36 4.47397 x10~* 4.19821 x107"” 2.60158 x107"” 2.58334 x107"” 

0.52 3.60599 x10~”? 3.48063 x10” 2.07758 x10? 2.06972 x10 

0.62 3.31023 x10 3.22403 x10~” 1.90174 x10°” 1.89643 x10°”? S 
0.81 2.96389 x10"? 2.91493 x10? 1.69762 x10~” 1.69467 x10~"? £ 


ARTICLE 


Table 8: Numerical values of Skin friction for various values of Prandtl number in H,O and C,H,.0O, based fluids when 


k, =1.380658 x10 JK |, € =1,0 = 25,d = 15nm, ¢ = 0.06mm* 


Pr tT (Al,O,)C,H,O, tT (Al,O,)H,O t(Cujye.4 0, tT (Cu)H,O 

0.71 —1,70014 x10” — 1.61226 x10" — 4.01659 x10” — 3.99391 x10” 
171 — 8.7099 x10"4 — 8.25961 x10"4 — 2.05783 x10'° — 2.04616 x10'° 
2.71  —5,47952 x10" —5.194.x10" — 1.29578 x10'4 — 1.28839 x10"4 
3.71 — 8.00715 x10” — 7.55025 x10" —1.91401 x10"° — 1.90237 x10"° 


Table 9: Numerical values of | skin friction for various values of nanoparticles volume fractions in H,O andC,H.O, based fluids 


when 

k, = 1.380658 x10 ~ JK ',€ =1,0 = 25,d =15nm,Pr = 0.71 

o(mm*) tT (Al,O;)C,H,O, tT (Al,O;)H,O tT (Cu)C,H,O, tT (Cu) H,O 
0.06 —1.70014 x10"’ — 1.61226 x10'’ — 4.01659 x10"” — 3.99391 x10" 
0.07 — 1.58285 x10"® — 1.49008 x10'® — 3.74798 x10"* — 3.72381 x10"® 
0.08 3.00003 x10"? 2.80444 x10!” 7.11892 x10'” 7.06775 x10'” 
0.09 2.0286 x10!” 1.88359 x10" 4.82372 x10" 4.78561 x10" 


Table 10: Numerical values of skin friction for various values of nanoparticles diameter in H,O andC,H.O, based fluids when 


k, =1.380658 x10 * JK |, € =1,0 = 25,6 =0.06mm’, Pr = 0.71 


d(nm) tT (Al,O,)C,H,O, tT (Al,O,)H,O tT (Cu)C,H,O, tT (Cu) H,O 

17 — 1.40912 x10" — 1.33629 x10" — 3.32914 x10" — 3.31026 x10" 
19 —1.19255 x10" —1.13095 x10" — 2.81757 x10" — 2.80159 x10" 
21 — 1.02634 x10" — 9.73293 x10'° — 2.4248 x10" — 2.41105 x10" 


Table 11: Numerical values of skin friction for various values of nanoparticles Sphericity in HO and C,H QO, based fluids when 
k, =1.380658 x10 * JK 1, é =1,0 = 25,0 = 0.06mm’ , Pr = 0.71, d(nm) = 15 
t (Al,O,)C,H,O, t (Al,O,;)H,O tT (Cu)C,H,O, tT (Cu) H,O 
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0.36 — 2.34957 x10" — 2.13572 x10" — 5.62845 x10" — 5.56934 x10" 


0.52 —1.70014 x10'’ — 1.61226 x10"’ — 4.01669 x10"’ — 3.99391 x10" 
0.62 —1.49532 x10'” —1.4373 x10'” — 3.51769 x10" — 3.50299 x10" 
0.81 —~1.2669 x10" —~ 1.23563 x10!” — 2.96682 x10'” — 2.95911 x10" 


Table 12: Numerical values of Thermal conductivity ratio for various values of nanoparticles volume fractions in H,O and 
C,H OQ, based fluids when 
k, =1.380658 x10 * JK |, € =1,0 = 25,d =15nm,Pr = 0.71 


: Nes kes ke. k 
mm 
ye) “ (Al,O;) —"_.(CxICjH.O, ” (cu) H,O ” (Al,O,)H,O 
f f Ky ky 
C,H 0, 
0.06 29.9320 29.9246 30,2434 30.2472 
0.07 12.6446 12.6358 -13.0116 -13,0162 
0.08 4.64418 4.65451 4.22025 4.21494 
0.09 21.9344 21.9462 21.4523 21.4462 


Table 13: Computed values of Thermal conductivity ratio for various values of nanoparticles Temperature in H,O and C,H .O, 


based fluids when 


k, = 1.380658 x10 ° JK ',€ =1,6 = 0.06mm’ ,d =15nm,Pr = 0.71 


0 Kn k k. ey 
—* cw C,H,0, a “ (cu) H,O —! (Al,O,)H,O 
Ky ky Ky Ky 

(Al,O,)C,H,O, 

25 29.9246 29.932 30.2434 30.2472 

50 -29,9246 -29,932 30.2434 30.2472 

100 -29,9246 -29,932 30.2434 30.2472 

1000 —-29.9244 -29.9318 30.2432 -30.2470 
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Table 14: Numerical values of Thermal Conductivity ratio for various values of nanoparticles diameter in H,O andC,H.O, 


based fluids when 


k, =1.380658 x10” JK |, € =1,0 = 25,6 =0.06mm’, Pr = 0.71 


d(nm) kup Kap kes 
— (Al,O,;)C,H,O, — (Cu)C,H,O, (Cu) H,O 
k k, k 

15 -29.932 -29.9246 -30.2434 

17 -29.932 -29.9246 -30.2434 

19 -29.932 -29.9246 -30.2434 

21 -29.932 -29.9246 -30.2434 


Table 15: Numerical values of Thermal conductivity ratio for various values of sphericity in C,H .O, and HO based fluids when 


k, =1.380658 x10 * JK ',é =1,0 = 25,d =15nm,¢ = 0,06mm’ Pr = 0.71 


¥ k,, j Kp k,, 
(Al,O;)C,H.O, (Cac, HO, (Al,O,)H,O 
k, k, k, 
0.36 -29.7771 -29.7618 -29.8246 
0.52 -29.9320 -29.9246 -29.9559 
0.62 -29.9890 -29.9837 -30.0061 
0.81 -30.0590 -30.0559 -30.0693 


Table 16: Numerical values of Temperature profile for various values of nanoparticles volume fractions in H,O andC,H .O, 


based fluids when 


k, =1.380658 x10 * JK ',é =1,0 = 25,d =15nm, Pr = 0.71 


d(mm*?) — @(r) (Al, 0,) C,H 05 @(r) (Al,O;)H,O A(r) (Cu) C,H,O, 
0.06 3.65787 x10 7! 3.65787 x10 7! 2.11613 x10~ 
0.07 9.29939 x10” 9.29939 x10” 5.26524 x10 
0.08 1.52039 x10” 1.52039 x10 6.88666 x10** 
0.09 4.20135 x10 7! 4.3694 x10! 2.41049 x10 7! 


Kp 
—(Al,O;)H,O 
k, 

-30.2472 

-30.2472 

-30.2472 


-30.2472 


Kap 
—(Cu)H,O 
k 

-29.7672 
-29.9272 


-29.9856 


-30.057 


Or) (Cu)H,O 


2.11513 x10! 


5.26524 x10 


9.18222 x10 °° 


2.53735 x10" 
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Table 17: Numerical values of Temperature profile for various values of Prandtl number in H,O and C,H .O, based fluids when 


k, = 1.380658 x10 ~ JK ',€ =1,0 = 25,d = 15nm, ¢ = 0.06mm* 


Pr Or) @(r)(Al,O,;)H,O Ar) (Cu)C,H,O, A(r) (Cu) H,O 
(Al,O,)C,H.O, 

0.71 3.65787 x10 ~~! 3.65987 x10 ~~! 2.11514 x10! 2.11513 x10! 

1.71 2.13685 x10” 2.9941 x10 1.21694 x10~°° 1.22402 x10~°° 

201 5.37676 x10 °*° 7.52176 x10 3.0611 x10°~° 3.0611x10°~° 

3.71 1.0055 x107'° 1.40866 x10” 5.74104 x10”? 5.74104 x10” 


Table 18: Numerical values of Temperature profile for various values of nanoparticles diameter in H,O andC,H,<O, based 


fluids when 

k, =1.380658 x10* JK |, € =1,0 = 25,d =15nm,Pr = 0.71 

a(nm) Ar) (Al,O;)C,H,O, @(r) (Al,O,)H,O O(r) (Cu)C,H,O, Or) (Cu)H,O 
15 3.65787 x10 ~~" 3.65787 x10! 2.11513 x10" 2.11513 x10" 
17 4.22552 x10! 4.14266 x10! 2.37682 x10~! 2.37682 x10 ~~" 
12 4.64235 x10 ~~" 4.64235 x10” 2.64499 x10”! 2.64499 x10 ~~" 
21 5.15683 x10! 5.15683 x10~! 2.91976 x10~! 2.64499 x10! 


Table 19: Numerical values of Temperature profile for various values of nanoparticles Sphericity in HO andC,H.O, based 


fluids when k, = 1.380658 x10-° JK", é =1,0 = 25,6 = 0.06mm’ , Pr = 0.71, d(nm) =15 


ue Or) (O(r) Al,O,) H,O O(r) (Cu)C,H,O,  O(r) (Cu) H,O 
(Al,O;)C,H.O, 

0.36 3.7357 x10 7! 3.7351 x10 7 2.11513 x10 2.11513 x10~ 

0.52 3.65787 x10 ~" 3.65787 x10 7! 2.11513 x10 2.11513 x10 

0.62 3.65787 x10 7" 3.65787 x10 7! 2.11513 x10* 2.11513 x10~ 

0.81 3.7357 x10 7! 3.65787 x10?! 2.11513 x10~ 2.11513 x10~ 
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Table 20: Numerical values of Velocity profile for various values of nanoparticles volume fractions in H,O and C,H ,O, based 


fluids when 
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k, =1.380658 x10 * JK |, € =1,0 = 25,d =15nm, Pr = 0.71, Re = 1000,G, = 0.92 


“p(mm?) u(r) (AL,0,)C,H,O, u(r) (AL,O,) HO u(r) (Cu) CHO, utr) (Cu) HO 
0.07 — 1.18714 x10"° —1.11756 x10"° — 2.81099 x10"° — 2.79286 x10"° 
0.08 2.25002 x10'! 2.10333 x10"! 5.33919 x10" 5.30081 x10'! 
0.09 1.52145 x10” 1.41268 x10” 3.61778 x10” 3.58919 x10" 


Table 21: Numerical values of Velocity profile for various values of Prandtl number in HO and C,H ,O, based fluids when 


k, = 1.380658 x10 JK, € =1,0 = 25,d =15nm,Re = 1000,G, =0.92 


Pr u(r) u(r) (Al,O,;)H,O u(r) (Cu)C,H,O, u(r) (Cu) H,O 
(Al,O,)C,H.O, 

0.71 —1.27511x10° —1.2092 x10° — 3.01253 x10° — 2.99544 x10° 

171 — 6.53311 x10° — 6.19543 x10° — 1.54349 x10’ — 1.53474 x10’ 

2.71 — 412363 — 391048 — 974235 — 968710 

3.71 — 62639.7 — 59401 .9 — 147991 —147151 


Table 22: Numerical values of Velocity profile for various values of nanoparticles diameter in H,O andC,H OQ, based fluids 


when 

k, =1.380658 x10 * JK |, € =1,0 = 25, Pr = 0.71, Re = 1000,G, =0.92 

d(nm) u(r) (Al,O,)C,H,O, u(r) (Al,O,)H,O u(r) (Cu)C.H0,. ue) (Cu)H,O0 
15 —1.27511x10° —1.2092 x10° — 3.01253 x10° — 2.99544 x10” 
17 — 1.05684 x10° —1.00222 x10° — 2.4968 x10° — 2.4827 x10° 
19 — 8.94446 x10° — 8.48214 x10° — 2.11319 x10” — 2.1012 x10° 
21 — 7.6976 x10° — 7.29973 x10° —1.81861 x10° —1.80829 x10° 


Table 23: Numerical values of Velocity profile for various values of nanoparticles Sphericity in H,O andC,H.O, based fluids 
when k, = 1.380658 x10 * JK |, € =1,0 = 25,d =15nm, Pr = 0.71, Re = 1000,G, = 0.92 
VY u(r) u(r) (Al,O;)H,O u(r) (Cu)C,H,O, u(r) (Cu) H,O 
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(Al,O;)C,H.O, 
0.36 — 1.76218 x10° 
0.52 —1.27511x10° 
0.62 —1.12151x10° 


—1.6018 x10° 


—1.2092 x10? 


— 1.07798 x10” 


— 9.26728 x10° 


— 4.22135 x10° 


— 3.01253 x10” 


— 2.63828 x10” 


= 999512 x10" 


— 4.17702 x10° 


— 2.99544 x10° 


— 2.62725 x10° 


0.81 — 9.50177 x10° 


— 2.21934 x10” 


Table 24: Numerical values of Velocity profile for various values of Reynolds’ number in H,O and C,H .O, based fluids when 


k, = 1.380658 x10 ° JK ',€ =1,0 = 25,d =15nm,Pr = 0.71,G, =0.92 


Re u(r) 


u(r) (Cu)C,H,O, 


(Al,0;)C, 8,0, 
1000 = =—],.27511x10° 


2000 =~ ].27511x10° 


3000 =9_].27511x10° 


u(r) (Al,O;)H,O 


— 1.2092 x10° 
— 1.2092 x10° 
— 1.2092 x10° 


— 1.2092 x10° 


— 3.01253 x10° 


— 3.01253 x10° 


— 3.01253 x10° 


— 3.01253 x10° 


u(r) (Cu) H,O 


— 2.99544 x10° 
— 2.99544 x10° 
— 2.99544 x10° 


— 2.99544 x10° 


4000 —1,27511x10° 


Table 25: Numerical values of Velocity profile for various values of Grashofs number in HO and C,H ,O, based fluids when 


k, =1.380658 x10” JK ',é =1,0 = 25,d =15nm,Pr = 0.71, Re = 1000, 


u(r) (Cu)C,H,O, 


G, u(r) 
(Al,O;)C,H.O, 
0.92 —1.27511x10° 
Le —1.27511x10° 
202 —1.27511x10° 


u(r) (Al,O;)H,O 


— 1.2092 x10° 
— 1.2092 x10° 
— 1.2092 x10° 


— 1.2092 x10° 


— 3.01253 x10° 


— 3.01253 x10° 


— 3.01253 x10° 


— 3.01253 x10° 


u(r) (Cu) H,O 


— 2.99544 x10° 
— 2.99544 x10° 
— 2.99544 x10° 


— 2.99544 x10° 


4000 —1,27511x10° 


Table 26: Computed values of Viscosity ratio of different shapes nanoparticles for various values of nanoparticles volume fractions 


Platelet 


g(mm*) 


0.06 5.43136 


Blade 


2.31988 


Cylinder 


5.06584 


Brick 
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2.81104 
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0.07 6.59874 2.62617 6.37656 3.44286 

0.08 7.88864 2.95712 7.86816 4.16896 

0.09 9.30106 3.31273 9.54064 4.98934 
4. DISCUSSION 


Mean Nusselt number 


Table 4 displayed the effect of Prandtl number on the heat transfer coefficient of Cu and Al,O, nanoparticles in H,O and 
C,H,O, based fluids. It shows that increase in Pr, corresponds to an increase in the nanofluids of Cu and Al,O, in both H,O 


and C,H .O,. However, the magnitude of the nanofluid of Cu is greater than nanofluid of Al,O, . Table 5, numerically showed 


that increase in nanoparticles volume fractions, enhanced the heat transfer coefficient of the nanoparticles in the two base fluids 
under consideration and this observation is consistent with the work of Polidori et al [26] and Abu-Nada [27]. The nanoparticles 
diameter is an important parameter in the enhancement of heat transfer coefficient and this fact is clearly demonstrated in Table 6. 
Almost all literature cited, established this fact. Increase in nanoparticles sphericity as depicted in Table 7, showed a decrease in the 
heat transfer coefficient. This is because mixture of nanoparticles with base fluids delayed diffusion and this in turn affect the Nusselt 
number. 


Skin friction 


Table 8, displayed that, as the Prandtl number increases, the skin friction within the nanofluids of Cu and Al,O; in H,O and 
C,H ,O, base fluids increases. However, Al,O, nanofluid increase in magnitude appreciates better than Cu nanofluid. 


Table 9, showed that the skin friction increases as the nanoparticle volume fraction increases but beyond 0.08mm? fraction, a 
decrease was observed. The explanation is that to improve the skin friction of nanofluids particularly Cu and Al,O, nanofluids, 
care must be taken to ensure the choice of the nanoparticle volume fractions. It is a truism that increase in skin friction between 
molecules nanofluids is enhanced by the size of the nanoparticles and Table 10 described it vividly in the two nanofluids under 
consideration. 

The skin friction increases as a result of increase in the sphericity of nanoparticles in the nanofluids. However, the increase is 
higher in water based fluids as shown in Table 11. This observation is consistent with the work of Temovera et al [28]. 


Thermal conductivity ratio 
Table 12 displayed the effect of increasing nanoparticle volume fractions on thermal conductivity of Cu and Al,O, nanofluids. It 
shows that increase in volume fractions led to an increase in both the water and ethylene glycol based fluids. This result is in 
agreement with that of Lee and Choi [29]. Table 13, showed that an increase in temperature up to about 373K does not show any 
appreciable enhancement in thermal conductivity of Cu and Al,O, nanofluids. This assertion was also corroborated by Li and 
Williams [30], Kolade et al [31] and Williams et al [32]. However, at temperature of 1273K and above, an increase in thermal 
conductivity was observed. 

An increase in nanoparticles diameter of Alumina and Copper in ethylene glycol and water based fluids as shown in Table 14, no 
increase in thermal conductivity was observed. From Table 15, an increase in sphericity results in a decrease in thermal conductivity 


enhancement but the decrease is more pronounced in C,H OQ, based fluid than HO based fluid, Lee and Choi [29] and Wang 
et al [20] also reported same findings. 

Temperature Profiles 

Table 16, showed numerical values which indicated that increasing the volume fractions of Cu and Al,O, nanofluids showed a 


decrease initially and later an increase in the temperature profile of the nanofluids but increasing the Prandtl number showed a 
corresponding increase in the temperature profile of the nanofluids as shown in Table 17. Table 18 showed that increasing the 
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diameter of nanoparticles also increases the temperature profile of the nanofluids. A decrease was observed in temperature profile 


of the nanofluids as the sphericity is increased as shown in Table 19. in the absence of the sphericity, the observations are 
consistent with the works of Koo and Kleinstrever ([18 and 19]), though weak owing to the presence of Brownian motion 


Velocity Profiles 
Tables 20 — 23, respectively displayed the numerical values of increasing volume fractions, Prandtl number, nanoparticles diameter 
and sphericity of the two nanofluids, showed that the velocity profiles of the nanoofluids all increased considerably but surprisingly 
an increase in Reynolds’ number(within the laminar and turbulence flow regime) and Grashof number due to temperature(which 
implies cooling the surface of the nanofluids) do not affect the velocity profiles of the nanofluids as shown in Table 24 and Table 25 
respectively. 

Table 26 clearly showed from the model that viscosity is largely dependent on nanoparticle volume fractions and therefore, its 
increase will no doubt corresponds to an increase in the viscosity of the nanofluids and it is supported by all the models listed in 
Wang and Mujumdar [33] and the work of Hashin and Shtrikman [34]. 


5. CONCLUSION 

The study theoretically tackled thermal conductivity enhancement and heat transfer coefficient of two nanofluids, namely alumina 
and copper in water and ethylene glycol based fluids. Although, the numerical values used showed considerable agreement with 
some existing works cited (Xuan and Li [35], Sato et al [36], a departure was also observed (Wong and Leon [37], Kuznetsov and 
Nield [38]) in others maybe as a result of approach or approximations. One vital observation is that a model however novel may not 
incorporate all the essential ingredients or factors of a given nanofluid. Temperature is also a weak parameter in influencing the 
enhancement of thermal conductivity. 


Funding 
This study has not received any external funding. 


Conflict of Interest 
The authors declare that there are no conflicts of interests. 


Peer-review 
External peer-review was done through double-blind method. 


Data and materials availability 
All data associated with this study are present in the paper. 


REFERENCE 

1. Choi S. U. S. Enhancing thermal conductivity of fluids with 6. Hamilton R. L, Crosser O. K: Thermal conductivity of 
nanoparticle, in: DA. Siginer, H.P. Wang (Eds.), heterogeneous two component systems. Industrial 
Developments and Applications of Non-Newtonian Flows. Engineering and Chemistry Fundamentals 1962, 1(3):187- 
ASME FED. (1995), 66:99-105 191. 

2. Yu, W. H, France, D. M, Routbort J.L and Choi S.U.S: Review 7. Jeffrey D. J: Conduction through a random suspension of 
and comparison of nanofluid thermal conductivity and heat spheres. Proceedings of Royal Society, A 1973, 335: 355-367. 
transfer enhancements. Heat Transfer Engineering 2009, 8. Davis R. H: The effective thermal conductivity of a composite 
29:432-460. material with spherical inclusions. International Journal of 

3. Eapen J, Rusconi R, Piazza R and Yip S: The classical nature of Thermophysics 1986, 7:609-620. 
thermal conduction in nanofluids. Journal of Heat Transfer 9. Wang, L. Q, Zhou, X. S, and Wei, X. H: Heat conduction 
2010, 132, 102402-1-102402-14. mathematical models and analytical solutions Berlin: 

4. Keblinski, P., Eastman, J. A., and Cahill, D. G. Nanofluids for Springer-Verlag; 2008. 
thermal transport. Materials Today, (2005). 8(6): 36-44 10. Koo J, Kang Y and Kleinstreuer C. A nonlinear effective 

5. Xie H, Wang J, Xi T and Liu Y: Thermal conductivity of thermal conductivity model for carbon nanotube and 
suspensions containing nanosized SiC particles. International nanofiber suspensions. Nanotechnology 2008, 19, 375705-1- 


Journal of Thermophysics. 2002, 23(2):571-580. 375705-7. 


PageDO 


ARTICLE 


11. 


12. 


13: 


14. 


15. 
16. 


ATs 
18. 
19. 
20. 


21. 
22. 


23. 


24. 


Jang S. P and Choi S. U.S: Effects of various parameters on 
nanofluid thermal conductivity. ASME Journal of Heat 
Transfer 2007, 129:617-623. 

Patel, H. E., Das, S. K., Sundararagan, T., Nair, A. S., Geoge, B., 
and Pradeep, T. Thermal conductivities of naked and 
monolayer protected metal nanoparticle based nanofluids: 
Manifestation of anomalous enhancement and chemical 
effects. Applied Physics Letters, (2003). 83, 2931-2933 

Das, S. K., Putta, N., Thiesen, P., and Roetzel, W. Temperature 
dependence of thermal conductivity enhancement for 
nanofluids. ASME Transnational Journal of Heat Transfer, 
(2003). 125, 567-574 

Xuan, Y., Li, Q., and Hu, W. Aggregation structure and 
thermal conductivity of nanofluids. AICHE Journal, (2003). 
49(4), 1038-1043 

Kumar, D. H., Patel, H. E., Kumar, V. R. R., Sundararajan, T., 
Pradeep, T., and Das, S. K. Model for heat conduction in 
nanofluids. Physical Review Letters, (2004). 93(14): 144,301- 
1-144,301-4 

Bhattacharya, P., Saha, S. K., Yadav, A., Phelan, P. E., and 
Prasher, R. S. Brownian dynamics simulation to determine 
the effective thermal conductivity of nanofluids. Journal of 
Applied Physics, (2004). 95(11): 6492-6494 

Putnam, S. A., Cahill, D. G., Braun, P. V., Ge, Z., and Shimmin, 
R. G. Thermal conductivity of nanoparticle suspensions. 
Journal of Applied Physics, (2006). 99( 8), 084,308 

Koo, J. and Kleinstreuer, C. A new thermal conductivity 
model for nanofluids. Journal of Nanoparticle Research, 
(2004), 6(6): 577-588 

Koo, J. and Kleinstreuer, C. Laminar nanofluid flow in micro- 
heat sinks. International Journal of Heat and Mass Transfer, 
(2005). 48(13): 2652-2661 

Wang, X. W, Xu, X. F and Choi, S. U. S: Thermal conductivity 
of nanoparticle-fluid mixture. Journal of Thermal physics and 
Heat Transfer 1999, 13:474-480. 

Keblinski P, Phillpot S. R, Choi S. U. S and Eastman J. A: 
Mechanisms of heat flow insuspensions of nanos-sized 
particles (nanofluids). International Journal of Heat and Mass 
Transfer 2002, 45:855-863. 

Tiwari, R. K and Das, M. K. Heat transfer augmentation in a 
two-sided lid-driven differentially heated square cavity 
utilizing nanofluids. International Journal of Heat and Mass 
Transfer. (2007). 50; 9-10. 

Asma, K; Khan, | and Sharidan, S. Exact solution for free 
convecton flow of nanofluids with ramped wall Temperature. 
The European Physical Journal-Plus. (2015), 130: 57-71. 
Aaiza, G; Khan | and Shafie S. Energy transfer in mixed 
convection mhd flow of nanofluid containing different 
shapes of nanoparticles in a channel filled with saturated 
porous medium. Nanoscle Research Letters. (2015). 10(490): 
1-14. 


25. 


26. 


27. 


28. 


29. 


30. 


eile 


32. 


33. 


34. 


35 


36. 


37, 


38. 


Mansour, R. B., Galanis, N., and Nguyen, C. T. Effect of 
uncertainties in physical properties on forced convection 
heat transfer with nanofluids. Applied Thermal Engineering, 
(2007). 27(1): 240-249 

Polidori, G, Fohanno, S, and Nguyen, C. T. A note on heat 
transfer modelling of Newtonian nanofluids in laminar free 
convection. International Journal of Thermal Sciences, 
(2007).46(8): 739-744 

Abu-Nada, E. Application of nanofluids for heat transfer 
enhancement of separated flows encountered in a backward 
facing step. International Journal of Heat and Fluid Flow, 
(2008). 29(1): 242-249 

Timofeeva, E. V; Jules, R. L and Dileep, S Particle Shape effect 
on Thermo Physical Properties of Alumina Nanofluids. 
Journal of Application Physics. (2009).106: 014304 

Lee S and Choi S. U. S: Measuring thermal conductivity of 
fluids containing oxide nanoparticles. Journal of Heat 
Transfer 1999, 121:280-289. 

Li, C. H, Williams W: Transient and steady-state experimental 


comparison study of effective thermal conductivity of 
AlL,O, -water nanofluids. Journal of Heat Transfer 2008, 


130, 042407-1-042407-7. 

Kolade, B, Goodson K. E and 
performance of nanofluids in a laminar thermally developing 
tube flow. Journal of Heat Transfer 2009, 13: 052402-1- 
052402-8. 

Williams, W, Buongiorno J and Hu, L. W: Experimental 


Eaton J. K: Convective 


investigation of turbulent convective heat transfer and 


pressure loss of alumina/water and zirconia/water 
nanoparticle colloids (nanofluids) in horizontal tubes. Journal 
of Heat Transfer 2008, 130: 042412-1-042412-7. 

Wang, X and Mujumdar, A. S. A review on nanofluids - part |: 
Theoretical and numerical investigations. Brazilian Journal of 
Chemical Engineering. 2008.25(4), 613 - 630 

Hashin Z and Shtrikman S: Conductivity of polycrystals. 


Physical Review 1963, 130:129-133. 


. Xuan, Y, and Li, Q. Heat transfer enhancement of nanofluids. 


International Journal of Heat and Fluid Flow. 2000, 21: 58-64. 
Sato, Y, Deutsch, E and Simonin, O. Direct numerical 
simulations of heat transfer solid particles suspended in 
homogeneous isotropic turbulence. International Journal of 
Heat and Fluid Flow. 1998. 19(2): 187-192 

Wong, K. V and De Leon, O. Application of Nanofluids: 
current and future. Advances in Mechanical Engineering. 
2000. 2: 519659. 

Kuznetsov, A. V and Nield, D. A. Natural 
boundary layer flow of a nanofluid past a vertical plate. 
International Journal of Thermal Science. 2010. 49: 243-247 


convective 


Pagel9 


